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ABSTRACT 

We study the inhomogeneous cosmological evolution of the Newtonian gravitational 
'constant' G in the framework of scalar-tensor theories. We investigate the differences 
that arise between the evolution of G in the background universes and in local inho- 
mogeneities that have separated out from the global expansion. Exact inhomogeneous 
sohitions are found which describe the effects of masses embedded in an expanding 
FRW Brans-Dicke universe. These are used to discuss possible spatial variations of 
G in different regions. We develop the technique of matching different scalar-tensor 
cosmologies of different spatial curvature at a boundary. This provides a model for 
the linear and non-linear evolution of spherical overdensities and inhomogeneities in 
G. This allows us to compare the evolution of G and G that occurs inside a coUapsing 
overdense cluster with that in the background universe. We develop a simple virialisa- 
tion criterion and apply the method to a realistic lambda-CDM cosmology containing 
spherical overdensities. Typically, far slower evolution of G will be found in the bound 
virialised cluster than in the cosmological background. We consider the behaviour that 
occurs in Brans-Dicke theory and in some other representative scalar-tensor theories. 
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1 INTRODUCTION 

Many of the most iriteresting extensions of the genera! the- 
ory of relativity are scalar-tensor theories of gravity. They in- 
corporate scalar fields which can carry space-time variations 
in scalar quantities that are traditionally assumed to be con- 
stant in general relativity. In this respect they provide nat- 
ural arenas in which to explore the consequences of varying 
'constants' of Nature, like Newton's 'constant' G, or Som- 
merfeld's fine structure 'constant', a. In this paper we shall 
focu s upon the first of these applications, pioneered by Jor- 
dan yordaril949l . ll955l . ll959l) and then refined into the most 
familiar generalisation of E instein's theory of gra vitation by 
Brans and Dicke in 1961, jBrans fc Dickelll96iri . This the- 
ory features both in direct explorations of the possible vari- 
ation of G, in dimensional redu ction of higher-dimensional 
theor ies, and in string theories JGreen. Schwartz fc Wittenl 
Il987l) . where is appears containing a dilaton with non- 
minimal coupling. The original theory of Brans-Dicke is 
recognised as the simplest case of a family of scalar-tensor 
gravity theories in which the original Brans-Dicke coupling 
parameter, tj, becomes a function of the scalar field that 
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carries variation s in G, JBergmannl 1196811 ■ iWagoneilll970l : 
lNordtveddll970lV 

Extensive studies have been made of co smological Solu- 
tions of scalar-tensor gravity theories llFuii fc Maeda 2003i) . 
although they are limited in two respects. First, they focus 
on the simplest case of isotropic expansion with zero spa- 
tial curvature, where simple exact solutions exist. Second, 
they are exclusively concerned with spatially homogeneous 
cosmologies. The latter restriction means that the value of 
G and its rate of change in time, G, are required to be the 
same everywhere in the universe. This assumption has run 
through the entire hterature on varying G and so it is gener- 
ally assumed, for example, that local observational bounds 
on varying G derived from geophysics, solar system dynam- 
ics, stellar evolution, or white dwarf cooling can be applied 
directly to constrain cosmological variations in G on extra- 
galactic scales or in the very early universe. There is no 
justification for this simp lifying assumption, as po inted out 
by Barrow and 0'Toole jBarrow fc O'TooldliÖoïh . The lo- 
cal bounds on varying G are all derived from observations 
of gravitationally bound 'lumps' which are in gravitational 
equilibrium and do not take part in the expansion of the uni- 
verse. Before they can be extrapolated to constrain possible 
variations of G in a background Friedmann universe (as is 
habitually done in the literature, without justification) we 
need to understand how G and G are expected to vary in 
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space in a realistic inhomogeneous universe. Since, even on 
the scale of a typical galaxy, the amphtude of visible den- 
sity inhomogeneities are of order 10 , we need to go beyond 
hnear perturbation theory for such an analysis. 

In this paper we begin to confront this deficiency by 
tracing the evolution of G, first in exact inhomogeneous Solu- 
tions and then a simple, but not unrealistic, inhomogeneous 
universe in which a zero-curvature Brans-Dicke-Friedmann 
background universe is populated by spherical overdensi- 
ties which are modelled by positive curvature Brans-Dicke- 
Friedmann universes in the dust-dominated era of the uni- 
verse's history. This will enable us to track the different evo- 
lution foliowed by G{t) in the background universe and in 
the overdense regions, which eventually separate off from the 
background universe and start to contract to high-density 
like separate closed universes. This process can produce sig- 
nificant differences between G and G in the background and 
in the overdensities. Eventually, the coUapse of the spheri- 
cal overdensities will be stopped by pressure and a compli- 
cated sequence of dissipative and relaxation processes will 
lead to viriahsation and some final state of gravitational 
equilibrium. This state will provide the gravitational envi- 
ronment out of which which stars and planetary systems 
hke our own will form, directly reflecting the local value of 
G{t) inherited from their virialised protogalaxy or its parent 
protocluster. The simple model we use for inhomogeneities 
in density and in G has many obvious limitations, notably 
in its neglect of pressure, deviations from spherical sym- 
metry, accretion, and interactions between inhomogeneities. 
Nonetheless, we expect that it will be indicative of the im- 
portance of taking spatial inhomogeneity into account in any 
attempts to use observational data to constrain cosmological 
models which permit varying G. It provides the first step in 
a clear path towards improved realism in the modelling of 
inhomogeneities that mirrors the route foliowed in Standard 
cosmological studies of galaxy formation with constant G. 

In section |21 we give the field equations and the field 
equations for scalar-tensor gravity theories. To fix ideas, we 
use some exact Brans-Dicke-Friedmann cosmological Solu- 
tions in section 13 to model the time variation of G in some 
new exact inhomogeneous solutions of the field equations 
which describe a spherical inhomogeneity in an expanding 
universe. In section|31we use exact solutions of flat and closed 
vacuüm Brans-Dicke-Friedmann cosmologies to illustrate 
the use of the spherical coUapse model for non-linear over- 
densities. We then apply the same techniques to dust-fiUed 
Brans-Dicke-Friedmann universes in Section |ïï| The numer- 
ical solution of the equations are presented and discussed in 
this section for Brans-Dicke and some other scalar-tensor 
theories. A summary of our principal results is given in|ü| 



2 SCALAR-TENSOR COSMOLOGIES 

The action for a scalar-tensor theory of gravity is given by 



üjidf) 



S=—- / d"a;V^(0ii+^f^(;^''9^<?!.9,<?i-H67r£„). (1) 



167r 

Here </) is a scalar field, R is the Ricci scalar, and Lm is the 
Lagrangian for matter fields in the space-time and the free 
function tj(<^) specifies the scalar-tensor theory. In Brans- 



Dicke theory, uj is constant. The action above is varied with 
respect to g^i, to give the field equations 



pM" „''"Pi ('„MP„''o- „A»i'_pa"\ j 



&1T 



rp)l 



(2) 



and with respect to to obtain the propagation equation 



00: 



2iü{(l>) + 3 



{8nT-Lü'{<j>)g^''<j>.aM, 



(3) 



where prime denotes differentiation with respect to 0. In this 
paper we will often consider Friedmann-Robertson- Walker 
(FRW) universes with the metric 



ds — dt — a (t) 



dr^ 



1 — kr^ 



+ r'^{de'^ +sm'^"''^^' 



where a{t) is the scale factor and k is the curvature param- 
eter. Using the FRW metric in eq. ^ gives 



2 é^ é 



Stt 



_ Stt (p — 3p) 
~ ~^ {2üj + 3) 



-3iï^- 



{2u + 3) 



and 



S-ïï 2 ,rr4> ^0 



(4) 



(5) 



(6) 



where H = a/a is the Hubble rate, over-dot denotes differ- 
entiation with respect to comoving proper time, t, p is the 
matter density, and p is the pressure. Each non-interacting 
fluid source p{p) separately satisfies a conservation equation: 

p + 3H{p + p)=0. (7) 

Substituting eqs. @ and ((üj into eq. Q gives 



H + H^ - H- + 



<j) Lü (j) 



Stt (3ptJ -I- 3p -|- pui) 



3<^ 



{2u + 3) 



-f 



2 {2u + 3) . 



(8) 



A general feature of the scalar-tensor field equations is 
that any solution of general relativity (hence u and <j) both 
constant) for which the energy-momentum tensor of mat- 
ter has vanishing tracé (eg vacuüm, black-body radiation, 
Yang-Mills field, or magnetic field) is a particular (<^ = con- 
stant) exact solution of the scalar-tensor gravity theory. A 
specification of cij(0) is required to determine the theory and 
close the system of equations. In general, we do not know 
the form of (^(iji) but if the theory is to approach general rel- 
ativity in an appropriate limit then we require both lü ^t oo 



and Lü'{(f>)uj 
limit. 



O to hold simultaneously in the weak-field 
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3 BRANS-DICKE COSMOLOGIES 



Consider the simpl est case of Brans-Dicke (BD) theory 
iJBrans fc Dickell961il to fix ideas about varying G. In these 
theories uj{4>) = lü is a, constant. The three essential field 
equations for the evolution of the BD scalar field (ji{t) and 
the expansion scale factor a{t) in a BD universe are ©, 
l|SJl, and 0. Now, lj is a constant parameter and the the- 
ory reduces to general relativity in the limit lu ^ oo where 
(p = G~^ — > constant. The form of the general solutions 
to the Friedmann metric in BD theories are fuUy under- 
stood (fBarrow'1992^, ('Gure vich. Finkelstein fc Rubanll973 
[Holden & Wands 1998). The vacuüm solution is the t — > O 
attractor for the perfect-fiuid solutions. The general perfect- 
fiuid solutions with equation of state 



p= (7"1)P 



(9) 



and k — O can all be found. At early times they approach the 
vacuüm solutions but at late time they approach particular 
power- law exact solutions |Nariai 1968); 



I^U\ ^ ^[2+2c^(2-7)]/[4+3^7(2-7)] 



Ht) = <t>ot 



[2(4-37)/[4+3l;7(2-7)1 



(10) 



(11) 



These particular exact power-law solutions for a{t) and 
4>{t) are 'Machian' in the sense that the cosmological evo- 
lution is driven by the matter content rather than by the 
kinetic energy of the free (j> field. The sign of (j) is determined 
by the sign of 4 — 87. 

These solutions are spatially homogeneous and so can- 
not teil us about the effects of any spatial inhomogeneity 
in cj> and p on observational tests of time- varying G — 4>~^- 
Next, we consider some simple exact inhomogeneous solu- 
tions of BD theory in order to gain some intuition about 
the likely effects of spatial inhomogeneity in G . We will 
find that these exact simple homogeneous solutions play an 
important role in determining the time dependence of G in 
inhomogeneous solutions. 



in a vacuüm under the transformation 
immediately obtain JFuü fc Maedall20(> 



r 



From (1121 we 



and 



R = 4>{R + 6Dr + öff^T.p F,, ) 



(13) 



where here R is the Ricci scalar, F = ln0~^" and DF = 
-^df_i{-\/—gg'^''dvV). To derive the Einstein field equations 
with a minimally coupled scalar field, we can extremise the 
Lagrangian density 



^ " ^^ ' 16^^ + g^"""^'^ '^"■ 



(14) 



to get Einstein's field equations, G^^ = SnT^,^, where 
T^u = V.M V".!- -èflMi'<?"''V'>a V'i/3- We now set </> = 
exp[ip / ^\ ] so that under the conformal transformation 
prescribed by II12II we find that (I14II becomes 



C^. 



aijTri^' R - enr - e^-'g^T,^ f, 



-IGvr 



1 / 3, 

+ Ï6^('^+2) 



~1 M" V^'M Vi" - 

9 -^^' 



(15) 



We see that DF disappears on integration by parts 
and so we can discard it. We also note that (/^"F,^ F,^ = 
é,fi 0,^, so that 11511 simplifies to 



iff"'^ 



'^m^^'^^'^ ^5'^''0.m'^.>')- 



(16) 



This is the Lagrangian density that gives the BD action Q, 

with Cmatter = 0. 

Starting with a solution of Einstein's field equations 
with a scalar field we can apply this conformal transfor- 
mation to arrive at a solution of the BD field equations in a 
vacuüm. A spherically symmetrie exact solution for the col- 
lapse of a minimally-coupled sc alar field, i/), in general rela- 
tivity is known and is given by JHusain. Martinez fc Nuned 
I199A 



3.1 An Inhomogeneous vacuüm Brans-Dicke 
Solution 

It is well known that BD theory is related to general rel- 
ativity through a c onformal transformation of the form 
JFuü fc Maedall2003h 



ds" = {qt + h){f\r)dt^ - r^(r)d^) 

-- R'^{r,t){de^ + sin 9d(j>'^), (17) 

where f{r) = (1 - ^^ , R\r,t) = {qt + b)r^(l - f )^-° 
and a = ±-^. The evolution of the minimally coupled scalar 
field, Vi in Ihe Einstein frame, is given by 



ff/j'^ ~ i9t^v 



(12) 



where (j) is the BD scalar field. Symbols with bars refer to 
quantities in the Einstein (general relativistic) conformal 
frame and symbols without bars refer to quantities in the 
Jordan (BD) conformal frame. This conformal equivalence 
allows us to exploit known solutions of Einstein's field equa- 
tions with a scalar field to find solutions to the BD field 
equations in a vacuüm. 

We proceed by first showing explicitly the conformal 
equivalence of general relativity with a scalar field and BD 



V'(r,i) = ± 



40F 



In 



d(l_|)^(,i + fe)v^ 



(18) 



Here, q, b, c and d are constants. Now under the transfor- 
mation 112II . where (f> — exp[ip / -^^], we obtain 



ds^ = _^iL! \A{rrdt^ - A{r)-"dr^] 

di//3A(r)«/^/3 '- -' 

^^^ {de +sin 



(19) 
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and 



,t)= [dA{r 



^'^B{t)^ 



1//3 



(20) 



where A{r) = 1-^, B{t) = qt + 6 and /3 = ±^/2lü + 3. Wo 
now assume that q ^ O (i.e. the metric is not static) and 

, 3_ %/3 

define the new time coordinate t — [qt + b)^ '^F . 

In the limit that c — > O the r-dependence of the met- 
ric is removed and the space becomes homogeneous. In this 
case we expect (I2Ü1 to reduce to the FRW Brans-Dicke 
metric given in the last section. We see from the form of 
12011 that the metric should reduce to that of a flat FRW 
Brans-Dicke universe. Insisting on this limit requires us to 
set (3 = V2tJ-f 3, q = j^^ and d 
metric: 



1. This leaves the 




1.1 



0.9 
0.8 



100 



G (r, t) 
G (2.1, 5) 



ds 



A{rT^^-vr,^d? 



-a(l+- 



-Jip-Vü) 



- A{r) -^f't 3/3-v^ X 

[dr^ + A{r)r^{de^ + sin^ 9d<j>'^)] . (21) 
Rewriting 1201 with these coordinates and constants gives 



lr,t)=(l-- 



72/(v^^~l) 



(22) 



A comparison of 1211 1 with 13211 shows that 1211 1 does 
indeed reduce to a flat vacuüm FRW metric in the limit 
c ^ O (an inhomogeneous universe requires c 7^ 0). The 
metric 12H is asymptotically flat and has singularities at 
t = O and r — 2c; the coordinates r and t therefore cover 
the ranges O ^ t < 00 and 2 ^ - < ex. 

The equations l|47^ and H22|l can now be used to con- 
struct a plot of G{r,t); this is done in FigureQwhich was 
constructed by choosing —jgi^^ 1221 1. From the form of 
G{r, t) we see that this choice corresponds to an overden- 
sity in the mass distribution (identified by comparison with 
the inhomogeneous Brans-Dicke solution with matter, found 
below). In this figure, uj was set equal to 100, and c set 
equal to 1. This plot shows how G can vary in space and 
time in an inhomogeneous universe which consists of a static 
Schwarzschild-like mass sitting at r = O in an expanding 
universe. As r — > O the solution approaches the behaviour 
of the static spherical vacuüm BD solution but as r ^ 00 it 
approaches the behaviour of a BD Friedmann universe. 



3.2 An Inhomogeneous Brans-Dicke Solution 

With Matter 

We now seek a solution of the Brans-Dicke field equations 
||5J with the form 

ds^ = e'dt^ - e^a^{dr'^ + r^dö^ + r^ sin^ öd^^) (23) 

where e" = e"*''', e^ = e''^''' and a = a{t). In Appendix 
^we show that a solution of the field equations ^ for a 
metric of the form 12311 is given by 



1 — 

2kr 
1 + -^ 



(24) 



Figure 1. This graph illustrates the possible space and time vari- 
ations that can arise in G in inhomogeneous solutions to the 
Brans-Dicke field equations, normalised at r = 2.1 and t = 5. 



(25) 



e^ = 


<^-^y 


/-, c X 2(fc-l)(fc+2)/fc 

G.ï) 




a{t) = ao 


2c.(2-^) + 2 
f t \ 3iJ7(2-7) + 4 



and 



{r,t) = <f)o 



for the matter distribution 





2(4-3t) 




t ^ 

to y 


, 3u>-,(2~-,) + i 


f^~2t 

U+2f7 



-2{k^-l)/k 



p{r,t) = po 



ao 



37 



1 — 

2kr 
1 -I- -^ 



(26) 



(27) 



(28) 



where k 



4,+2bj 
3+2lü ■ 



This separable solution displays the same 



time dependence as the power-law FRW Brans-Dicke uni- 
verses, l)10^ - l|ll|l . but with an additional inhomogeneous r- 
dependence created by the matter source at r = 0. Such a 
distribution of matter in space is illustrated by figure|21 Here 
we have chosen, for illustrative purposes, 'y = 1, uj = 100 
and a background value set by the choice p — Pfrw (pfrw 
being the matter density that would be expected in the cor- 
responding homogeneous universe) . The temporal evolution 
of p is exactly the same as the FRW case. 

We see from Figure|5|that the matter density is isotropic 
and asymptotically constant as r ^ 00 with a sharp power- 
law peak near the origin. Now 1271 gives us 



G(r,i) = Go 



1 ol " 


r 


-!)/* 2(4-3^) 


1+2^7. 


l 3u,-,(2--,) + 4, 



(29) 



Equation 1291 is used, with the values ui = 100, 7 = 1 
and c = 0.5 to create Figure |21 which shows the space-time 
evolution of G{r,t). 

These results show how G{r, t) can vary in space and 
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Figure 2. Distribution of p as a function of r, with w = 100, 
from eq. fg<SD and c = 0.5. 



more slowly than the background, before reaching an expan- 
sion maximum and coUapsing back to high density, whilst 
the background continues to expand. In this section we con- 
sider vacuüm uriiverses only, so there exists spherically sym- 
metrie inhomogeneity in the expansion rate and in <ji ~ G^^ , 
but p = p — 0. 

For flat vacuüm FRW universes eq. Q gives 



a\ a(t)b_ _ i£ f (j>b_ 
aj a <f)b Q \4>b 



(30) 



where ' = 4- and 

at 



b = (pb(t) is the BD scalar field and a is the 
scale factor in the flat background. For a positively curved 
(fe = +1) region the scale factor is taken to be S{t), which 
satisfies the Friedmann equation for the closed vacuüm BD 
universe: 



G (r, t) 
G (1, 5) 




5" 



S'c^'p 



S I S ó„ 6 



g2' 



(31) 



where ' — ■£:, (j>p = (/f>p(T") and r and k are the proper time 
and curvature of this perturbed region. 

In matching these two regions at i = to = to we must 
satisfy the boundary conditions 



S{to) = a{to), 
(ppiTo) = 06 (io) 



and 



dS_ 
dr 

d(j>p 

"d7 



da 
H 

d4b_ 
dt 



100 



Figure 3. Evolution ofG{r,t) in space and time in an inhomo- 
geneous matter dominated Universe, with w = 100, from eq. 1296 
and c = 0.5. 



time in an asymptotically-flat universe with a peak of mat- 
ter at the origin. Observers located near the mass concen- 
tration will determine different values of G locally although 
they will find the same values of G/G everywhere because of 
the separable nature of the G{r, i) evolution in eq. 1281 . This 
was also the case for solution 11221 given in subsection 13.21 
In the next section we shall consider a more realistic model 
in which both G and G/G are different from place to place. 
Plots like Figure|21can be generaled for universes dominated 
by other types of cosmological fluid and with different rates 
of density fall off with r. 



4 MATCHING TWO VACUÜM FRW 
BRANS-DICKE UNIVERSES 

We will now consider a simple model of a spherically sym- 
metrie cosmological inhomogeneity produced by matching 
together flat and positively curved vacuüm FRW-BD uni- 
verses. This is a well studied technique, first introduced by 
Lemaitre, for studying the non-linear evolution of overden- 
sities in general relativistic FRW universes. The overdense 
region is modelled as a closed universe that at first expands 



4.1 The Background Universe 

Assuming solutions of the form (f>b oc t^ and a (x t^ and 
setting a(0) = O gives, on substitution into 13011 and 0, the 
fc = O BD vacuüm solutions (0'Hanlon fc Tupper 19721) 



^U\ ^ ^^ (1+2(1- y/3{3+2üj)r^) 



(32) 



and 



(l>b{t) = 't'boi-r) 



N -2(1-^3(3+2lj))-^ 



4.2 A CoUapsing Universe 

For the closed region we now foUow the m ethod given in refs. 
iBarrowlll992l : iBarrow fc Parsonall997r) to find expressions 
for S{t) and 4>v{t)- We start by introducing conformal time, 
77, deflned by Sdrj = dr; then eq. Q becomes 



t>p,r)r) + „"5,7; <Pp,7ï — U. 



This integrates direclly to yield 



V2,A{2üj + ?,)-^'^ 



(33) 
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where ^ is a constant. We now introducé the variable y 
(f>pS^ to write ((üj as 



R(t), 



2 



-4ky^ + -(bp,lS\2u: + 3). 



(34) 



Now, eqs. 1331 and 1341 give 



^ = VSAy-\2Lü + 3)-^ 



/2 



and 



y,l=-4ky^+A\ (35) 

The Solutions of eqs. 11351 . when k > O, are given by 

yiv) = ^sm{2Vk{v + B), (36) 



and 



Mv) = Ctanv/ (2-+3) (^^(^ + B)) (37) 

where B and C are arbitrary constants. We now fix the 
conformal time origin by setting i? = O, so that y — (j>pS'^ 
gives 



S{v)(x 



.1/2 



{2Vkv) 



tanvZ 4(2^+3) (^^) 



(38) 



4.3 From ri to t 



The function T{rj) is now obtained by integrating Sdrj = dr 
and T{rj) can then be used to obtain S{t). We now re- 
quire a relation between t and r, for this we proceed as 
in ref. JBarrow fc Kunzd ll^^7^ an d use the equation o f 
relat ivistic hydrostatic equilibrium (IHarrisonI ll97(J . Il973l) . 
JLandau fc Lifshit z 1975) 



0$ 
dr 



dp/dr 
P + P 



(39) 



where $ is the Newtonian gravitational potential and r is 
radial distance. Eqn. H39^ is derived under the assumptions 
that the configuration is static and the gravitational field is 
weak, so $ completely determines the metric; then 1391 is 
given by the conservation of energy-momentum for a perfect 
fluid. We now use dr = e dt, and for a scalar field we have 



an effective state with p = p and p - 
results gives 



dr 
H 



4b (t>p 

è' ^1/2 ' 



Combining these 



(40) 



Now 06 ex a '' and (f>'p cc S '' , and so with 1371 and 13211 
this gives 



dr _ _sin^/2(2\/fcr;) S^ a^^+i^ 
dt ^ 



sin 



'^^{2Vkrio) Si ^VïTF-4 



(41) 



-'o 



Eq. 1411 and the relation Sdrj — dr allow us to obtain 




Figure 4. Evolution of the scale factor S in the perturbed over- 
dense region (dashed line) and in the background afsolid line) 
with respect to the comoving proper time in the flat background. 




Figure 5. Evolution of G{t) in the overdense perturbed overdense 
region of positive curvature (dashed line) and in the spatially flat 
background universe (solid line). 



S{t) from H38|l . This is done numerically. Now fixing the 
constants of proportionality together with k in eqs. 1381 . 
1371 and 1321 . in order to satisfy the boundary conditions, 
we find equations for the evolution of the scale factors and 
scalar fields in the fiat background and perturbed region. 
These are matched at a boundary, at time to = tq = tjq. 



4:A Results 

Figure |1] shows the evolution of a{t) and S{t) when the re- 
gion described by S{t) becomes positively curved at initial 
time to = 1. In Figure |1| we choose ui = 100, for illustrative 
purposes, and ao = 5*0 = 1 so that the boundary condition 
for the matching of the first derivatives of the scale factors 



is given by 



> dt 1 



We can now express G — G(t) in the regions of dif- 
ferent curvature using equations 1471 . 1371 and 1321 . along 
with the appropriate coordinate transformations. This gives 
Figure 121 below. It is clearly seen that the evolution of G{t) 
is quite different in the two regions, as expected. The col- 
lapsing overdensity evolves faster than the background and 
possesses a smaller value of G but a larger value of \G/G\ at 
all times after the evolution commences. 
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5 A MORE REFINED SPHERICAL COLLAPSE 
MODEL 

We now generalise the spherical coUapse model described 
in the last section to the more astronomically realistic case 
of a flat universe co ntaininR matter and a cosmoloeical c on- 
stant (see. e.R. refs . JPadmanabhanlf ^95^ . (|PeacoclJl999() or 



while ^ reduces to 



JLahav et allf993) 'l. As before, we match a flat Brans-Dicke 
FRW background to a spherically symmetrie overdensity at 
an appropriate boundary and allow the two regions to evolve 
separately. 



5.1 The background universe 

Again, we consider a flat (fc = 0), homogeneous and isotropic 
background universe. Since we are interested in the matter- 
dominated epoch, when structure formation starts, we can 
assume that our universe contains only matter and a vacuüm 
energy contribution so that p — pm + Pa and p = pa = —pA 
give the total density and pressure, respectively. So, for the 
flat background, eq. gives for a general iü{(j>) theory. 



+ 



a (j) 6 



uj(j>^ Stv, , 



and eq. @ gives 



4> + 3U=, ^" . 



{pm + 4pa) - 



UJ(j} 



{2uj + 3) 



(42) 



(43) 



Here, pm oc a{t)^ and pA = constant. These equations 
govern the evolution of (l>{t) and a{t) in the flat expanding 
cosmological background. The contribution of the vacuüm 
energy stress (p = —p) to the Friedmann equation II42I I in 
BD cosmology differs from that in general relativity because 
of the presence of the variable field: it is not the same as 
the addition of a cosmological constant term to the Rif S of 
14211 . lïowever, with this proviso, we shall continue to refer 
to Acdm models in Brans-Dicke theories in the following 
sections. 



5.2 The overdensity 

Again, we consider a spherical overdense region of radius S 
and model the interior space-time as a closed FRW Brans- 
Dicke theory universe, ignoring any anisotropic effects of 
gravitational instability or coUapse. As usual, we assume 
there is no shell-crossing; this implies mass conservation in- 
side the overdensity and i ndependence of the radial coordi- 
nate iPadmanabharll995h . The evolution equations can now 
be written in a form that ignores the spatial dependence of 
the fields. Put p = pcdm + Pa and p = pa = — pA, where 
'cdm' corresponds to cold dark matter, so that (|HJ gives 



S-S^ ^ 



-^It^ 



2 (2lüc + 3) (j}c 
Stt (pcdm(3 -I- ujc) + Pa{3- 2lüc)) 



3(j>c 



{2üJc + 3) 



-^■^S"^'' (2a;, + 3) 



{pcdm + 4Pa) 



(45) 



(44) 



i2ujc + 3) 

where S = S{t) is the scale factor and 4>c ~ 0c (t) is the 
BD scalar field in the coUapsing region of positive curvature 
where Pcdm oc S{t)~^ and pA = constant. These equations 
give the evolution of 0c (t) and S{t). 

We have assumed that the equation of motion of the 
field inside the cluster overdensity is described by the local 
space-time geometry. This means that the field foUows the 
dark-matter coUapse from the beginning of the cluster's for- 
mation. We do not consider this to be fuUy realistic since 
there is expected to be an outfiow of energy associated with 
(j) from the overdensity to the backgr ound universe, as first 
notice d by Mota and van de Bruck JMota fc van de Bruckl 
||200^). The details of this outfiow of energy and its effect 
on the collapse can only be determined by a fuUy relativis- 
tic hydrodynamical calculation, which is beyond the scope of 
this study. Nevertheless, at late times during the collapse of 
the dark matter (and especially when the density contrast in 
the dark matter is very large) the field should no longer feel 
the effects of the expanding background and will decouple 
from it. We are also neglecting the effects of deviations from 
spherical symmetry, which grow during the collapse in the 
absence of pressure, along with rotation, gravitational tidal 
interactions between different overdensities, and all forms of 
non-linear hydrodynamical complexity. 



5.3 Evolution of the overdensity 

Consider the spherical perturbation in the cdm ffuid with a 
spatially constant internal density. Initially, this perturba- 
tion is assumed to have a density amplitude Si > O where 
\Si\ ^ f. The initial cdm density inside the overdensity is 
therefore pcdm = Pm{l + Si). 

Four characteristic phases of the overdensity's evolution 
can be identified: 

• Expansion: we employ the initial boundary condition 
(f)c = (j) fl'iid assume that at early times the overdensity ex- 
pands along with the background. 

• Turnaround: for a sufficiently large St, gravity prevents 
the overdensity from expanding forever; the spherical over- 
density breaks away from the general expansion and reaches 
a maximum radius. Turnaround is defined as the time when 
S = Sma^,S^Q!mAS <Q. 

• Collapse: the overdensity subsequently coUapses (5* < 
0). If pressure and dissipative physics are ignored the over- 
density would collapse to a singularity where the density 
of matter would tend to infinity. In reality this singular- 
ity does not occur; instead, the kinetic energy of collapse is 
transformed into random motions. 

• Vtrialisation: dynamical equilibrium is reached and the 
system becomes stationary with a fixed radius and constant 
energy density. 

We require our spherical overdensity to evolve from the 
linear perturbation regime at high redshift until it becomes 
non-linear, coUapses, and virialises. Thereafter, the overden- 
sity will become gravitationally stable and further local evo- 
lution of the scale factor and of the scalar field will cease. 
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However, the background scale factor will continue to ex- 
pand, and SC the background density and background value 
of G wiU continue to decrease. As a result, a significant dis- 
parity between the evolution of G inside and outside the 
cluster can result. 



5.4 Virialisation 

In scalar-tensor theories we expect that the gravitational 
potential will not be of the Standard local r~^ form. This 
requires reconsideration of the virial condition. According 
to the virial th eorem, equilibrium will be reached when 
jGoldsteinllïÖsÖI) 






(46) 



where T is the average total kinetic energy, U is the average 
total potential energy and S here denotes the radius of the 
spherical overdensity. 

The potential energy for a given component x can 
be calculated from its g eneral form in a spherical region 
jLandau fc LifshitjIÏÖTa) 



Ux =2% ptot^xT dr, 
Jo 

where ptot is the total energy density and ^^ is the gravita- 
tional potential due to the density component p^:- 

The gravitational potential ^x can be obtained from 
the weak-field limit of the field equations ||5J . This results in 
a Poisson equation where the terms associated to the scalar 
field can be absorbed into the definition of the Newtonian 
constant as ^Willlll993^ 



G^ = 



■i + 2cj,{(l)c) 1 



(47) 



3 + 2lüc(4>c) 4>c 
This results in the usual form for the Newtonian poten 



$,(s) = -27rG,p,(37,-2) S' 



tial 



where Ge is given by equation 14711 and 7^, — 1 is Px/ px for the 
fluid component with density px and pressure px (appear- 
ing due to the relativistic correction to Poisson's equation: 
A'ï> = 47rG(p + 3p)). 

In Acdm models of structure formation it is entirely 
plausible to set 7^, = 1 as the energy density of the cosmolog- 
ical const ant is negligible on the viri a lised scales we are con- 
sidering JWang fc Steinhardd Il998l : iMota fc van de Bruclj 
|2004). The potential energies associated with a given com- 
ponent (a;) inside the overdensity are now given by 



Ux 



., r. ^CptotpxO . 
15 



(48) 



Therefore, the virial theorem will be satisfied when 






where 

dU_ 
dS 



G, 



dS 

ve V 

PA + P4,a , P0C = "^^ 



dptot q5 , ^ ^ dpx c.5 I ^ „ „ c; c4 
PxO + GcPtot-^TTJ-O + GcptotPxOb 



dS 



(49) 



and we have used Utot = Ucd,n + Ua + U^^ , ptot = Pcdm + 
h^/è and 



dGc 
dS 



Ge 

S 



-Ge 



3 + 2ljc 
'4 + 2^^ 



Gc + 



2lü, 



c) 



(3 -I- 2a;, 



The other components of eq. 11491 are obtained from eqs. 
int and iBïïl . 

Using eq. II46L together with energy conservation at 
turnaround and virialisation, we obtain an equilibrium con- 
dition in terms of potential energies only 



1 (du 

2^""-U5 



+ Utot{Zv) = Utot{zta), 



(50) 



where z^ is the redshift at virialisation and Zta is the redshift 
of the over-density at its turnaround radius. The behaviour 
of G during the evolution of an overdensity can now be ob- 
tained by numerically evolving the background eq. 142t - 143j 
and the overdensity eqs. (I44II and 1451 until the virial con- 
dition 1501 holds. 

We point out here an inconsistency when one makes use 
of equation 1501 together with the assumption that energy is 
not conserved. This inconsistency is removed by assuming a 
negligible outflow of from the overdensity, in which case we 
regain energy conservation within the system and so retain 
self-consistency. 



5.5 Overdensities vs Background 

5.5.1 Brans-Dicke theory 

The Brans-Dicke coupling parameter ui is constant an d con- 
strained by a variety of local gravitational tests (see ijUzarJ 
'2003^ for a review). The strongest constraint to date is de- 
rived from observations of the Shapiro time delay of sig- 
nals from the Cassini space craft as it passes behind the 
S un. Th ese considcrations led Bertotti, less and Tortora 
(Bertptti. Icss & Tortora 2003), after a complicated data 
analysis process, to claim that u must have a value greater 
than 40000 (to 2a). This limit on uj must be satisfied at all 
times in all parts of the universe, and leads to the conclusion 
that Brans-Dicke theory must be phenomenologically very 
similar to general relativity throughout most of the history 
of the universe. However, we do still expect a cosmologi- 
cal evolution of the Brans-Dicke field (j> which determines 
the value of Newton's G; and we expect this evolution to 
be different in regions that coUapse to form the structure 
probed by Cassini compared to that in the idealised expand- 
ing cosmological background, as described above. Hence we 
expect the measurable value of G to be different in these 
two distinct regions with different histories. It is quite pos- 
sible that the region of gravitational equilibrium probed by 
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Figure 6. Plots of G against ln{l + z) in the background (dashed-line) and in an overdensity (solid-line), for different values of ui. 
Initial conditions are chosen in both cases so as to give G = Go, the present value of the Newton constant, at virialisation. We note 
that increasing u> decreases the difference in G between the overdensity and the background. 



Cassini and other local observations would find no percep- 
tible change in G locally despite the presence of change on 
cosmological scales outside of bound inhomogeneities. 

In this section we quantify this difference in G by nu- 
merically evolving a, S, (f> and (fic until virialisation occurs. 
At virialisation, the evolution of 5* and (f)c is expected to 
end, giving a value of G that is locally constant in time even 
though the cosmological evolution of a and continues. 

The plots in figure l^were constructed using the repre- 
sentative values ui = 40000 and lj = 500 and the boundary 
condition ^co — Gq^, so that the value of G measured inside 
the overdensity at present is equal to the value of Newton's 
G, as measured locally. The evolution of the background was 
determined by matching 0c to (f> at the time when the over- 
density decouples from the background, ti. We see a clear 
difference in the evolution of G in the two regions, as ex- 
pected. This example shows that we expect different values 
of G and G/G inside and outside virialised overdensities. 
The present value of G and G/G depends on the history of 
the region where it was sampled, as well as on the Brans- 
Dicke coupling parameter, lü. 

It can be seen from the plots in figure |ü] that increasing 
Lü has the effect of decreasing the difference in G between 
the background universe and the overdensity. The size of 
this inhomogeneity is found to be of order 1/lü and, corre- 
spondingly, reduces to zero a& uj ~* oo. This is an important 
consistency check for the methods used as we expect Brans- 
Dicke theory to reduce to general relativity, with a constant 
G, in this limit. 



5.5.2 Scalar-tensor theory with 2(j; + 3 



2A 1-^ 



choice of coupling function 



2üj{<f)) + 3 = 2A 



4>a 



Next, we consider a scalar-tensor theory with a variable 
üj{(j)). We investigate the class of theories defined by the 



where A, (j>ao, and p are positive definite constants. We re- 
fer to this as Theory 1. Such a choice of coupling was con- 
sidered by Barrow and Parsons and was solved exactly for 
the case of a flat FRW u niverse containing a perfect fluid 
jBarrow fc Parsonslll997fl . 

Setting the constants as 2A = {4>oo/ P) and p — 2 gives 
us th e scalar-tensor theory con sidered by Damour and Pi- 
chon JDamour fc Pichon ||l99 9|) and by Santiago, K alligas 
and Wagoner iSantiaeo. Kalligas fc Waeoneil [1993) • This 
choice of i.ij{(j)) corresponds to setting In A{<l)) — lnA{(j}o) + 
|/3(0oo — 4')^y where A^{4>) is th e conformal factor l/é from 
eq. H12^ . Damour and Nordvedt iJDamour fc Nordvedtll99.'tl 
consider this function as a potential and therefore justify 
its choice in relation to the generic parabolic form near a 
potential minimum. Expecting the function to be close to 
zero (i.e. GR), Santiago, Kalligas and Wagoner justify its 
expression as a perturbative expansion. This choice oi uj{(j)) 
with p > 1/2 corresponds to a general two-parameter class 
of scalar-tensor theories that are close to GR and will be 
drawn ever closer to it with cu -^ oo and lü' /uj^ ^ O as 
the universe expands and -^ (f>aa- We therefore consider 
it as a representative example of a wide family of plausible 
varying-G theories that generalise Brans-Dicke. 

The evolution of this form of 0^(0) is shown graphically 
in figureQfor difi'erent values of A and p. Clearly the evolu- 
tion of oj{(f)) is sensitive to both A and p and so the choice of 
these parameters is important for the form of the underlying 
theory. For illustrative purposes we choose here the values 
p = 1.5, 2 and 5 and A = 1,2 and 5. 

In a similar way to the BD case we now create an 
evolution of ^c that virialises at z = O to give the value 
Gco = 6.673 X 10~^^ , as observed experimentally. The corre- 
sponding evolution for 06 is calculated as before by match- 
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(a) 2Lü{(f>) + 3 for p=2 and A=l, 2 and 5 (solid, dashed 
and dotted Unes, respectively) . 



(b) 2uj((j>) + 3 for A=l and p=1.5, 2 and 2.5 (solid, dashed 
and dotted lines, respectively). 



Figure 7. In these gravity theories there is fast approach to general relativity at late times when 
behaviour at early times. 



but significantly different 



ing it to the value of (pc at the time the overdensity decou- 
ples from the background and begins to coUapse. In creating 
these plots we have used the conservative parameter values 
p — 2, ujco ~ 1.2 X 10^ and A = 6x 10~^ which are consistent 
with observation and allow structure formation to occur in 
a similar way to general relativity. The results of this are 
plotted in figure |H| 

Again, we note the different evolution of G{t) in the two 
regions, and the difference in the asymptotic values of G. We 
note that experimental measurements of G on Earth have a 
significant uncertainty with the 1998 CODATA value carry- 
ing an uncertainty 12 times greater than th e Standard value 
adopted in 1987. Th e 1998 value is given as llMohr fc Tavloj 
l200d:IScherTë3ll999^ 

Gi998 = 6.673 ±0.010 X 10"** cm^ffm~^s"^ 

while the 2002 CODATA pre-publication announcement 
reverts to the earlier higher accuracy consensus with 
JNational Institute of Standards and Technologvll2002f) 



G2002 = 6.6742 ± 0.0010 x 10 



cm gm 



We could re-run the above analysis with different values 
of A and p, but expect that the results would look qualita- 
tively similar. From Figure |7] we see that increasing (de- 
creasing) the values of A and p will increase (decrease) the 
value of ü;((/)) for a given (j>, thereby making the theory more 
(less) like GR. We therefore expect an analysis with a higher 
(lower) value of A and/or p to look very similar to the anal- 
ysis presented above with a less (more) rapid evolution of 
G{t). For the sake of brevity we omit such an analysis here. 



5.6 Space and Time variations of G 

We now calculate how time and space variations of G evolve 
with redshift and depend on the cdm density contrast, Ac. 
In order to do this, we make the definitions 



AG, , 



G{t) - Go 
Go 



and 



f« 



i(Zt.) 



Gc(t)-Gft(t) 
Gb{t) 



Pm \^v ) 

where Ge and Gt correspond to G as measured in the over- 
density and in the background universe respectively. 

The results of our numerical calculations, for a cluster 
which virialises at z„ = O, are presented in Figures |ÏÏ| 1101 
and 1111 respectively. These plots display the evolution of 
G/G, AG/G, and SG/G with redshift for Brans-Dicke, the 
theory of subsection 15.5.21 and some other choices of ui{(j)) 
that are specified in the captions. The parameters used in 
generating these plots are B = 0.4, G = 10~^®, D = 80, 
y4 = 6 X 10~^, p = 2 and lü = A x 10*^. These values were 
chosen so as to agree with observation and so that struc- 
ture formation is not significantly different from that which 
occurs in general relativity. 

It is clear from the plots that different scalar-tensor the- 
ories lead to different variations of G. The predictions of 
these models can be quite diverse. While some models pro- 
duce higher values of G inside the overdensity, others pro- 
duce a lower one. A feature common to all models is that 
the value of G and G/G inside an overdensity is different 
from G and G/G in the background universe. The reason 
for these differences is that in the non-linear regime, when 
the overdensity decouples from the background expansion at 
turnaround, the field <j) that drives variations in the Newto- 
nian gravitational "constant" stops feeling the background 
expansion. After turnaround, the field inside the overdensity, 
(jjc, deviates from the field, 0, in the background universe, 
leading to spatial variations in G. In reality, such spatial in- 
homogeneities in the value of G are small: SG/G ~ 10~^, 
Figure 1111 The time variations of G are even smaller than 
the spatial inhomogeneities but with a marked difference 
between the inside and outside rates of change. We find 
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A graph of the variation G against ln{l + z) for the background universe (dashed-line) and inside 



Figure 8. 2c<; + 3 = 2A 1 - -r^ 

an overdensity (solid-line) which give G = Go, the presently observed terrestrial value, at virialisation. 




log(1+z) 



Figure 9. Evolution of AG/G as a function of logil + z) for 
overdensities which virialise at z = O in a Acdm model. Upper 
plot: evolution inside the overdensity. Lower plot: evolution in the 
background universe. Thick solid line 2c<; + 3 = 4 x 10^, thin solid 
line 2a;+3 = B'^(p, dashed-line 2üj+3 = 2A\l--^\-^, dash-dotted 
-4 j„ii„j ;.'„_ o,., I q _ ^t\^ _ /• \2|-1 



line 2w + 3 = C|ln{^)|-*, dotted line 2(^ + 3 = D|l 

. Each model is normalised in order to have Gq = Gc{z 

inside the overdensities. 
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Figure 10. Evolution of G/G as a function of log{l + z) for 
overdensities which virialise at z = O in a Acdm model. Upper 
plot: evolution inside the overdensity. Lower plot: evolution in the 
background universe. Thick solid line 21^; + 3 = 4 X 10^, thin solid 
line 2a;+3 = B^ó, dashed-line 2üj+3 = 2A\l — -f-\~^ , dash-dotted 



2A|1-^| 



line 2c<;-f 3 = C|ln{^)|-*, dotted line 2lj -f 3 = D\l 
. Each model is normalised in order to have Go = 
inside the overdensities. 



- ( 't' '|2|-l 

Gc{z = 0) 



Gc/Gc ^ 10"^°s"^ inside the clusters and G/G sC 10""s"^ 
in the background, Figure HÜl 

Figures 1121 1131 and 1141 represent the values of AG/G, 
5G/G and G/G at the redshift of virialisation (zv). Once 
again the differences between the different scalar-tensor the- 
ories and between Gc(zv) and G{zv) are evident. 

Figure 1151 shows how the cdm contrast, Ac = 
pcdm{zv) / Pm{zv) , affccts the difference between the value 
of G inside an overdensity and in the cosmological back- 
ground. (Recall that, in the Einstein-de Sitter model 
Ac = Pcdm{zv) / Pm{zv) ~ 147 at virialisation, and Ac = 
Pcdm(zv) / Pm{zc) ~ 187 at coUapse) . It is interesting to see 



that different scalar-tensor theories produce a different de- 
pendence. 

From the numerical simulations, we see that a varia- 
tion of the Newtonian gravitational constant, of the order 
presented here, does not affect the predictions of the struc- 
ture formation models. For instance, the virialisation radius 
and the density contrast Ac of the virialised clusters are 
very similar to the Acdm and Standard cdm models with 
constant G . Nevertheless, one should point out that differ- 
ent scalar-tensor theories lead to slightly different results. 
Besides that, there is also a dependence on the initial con- 
diti ons, as with many other cosmological scalar fields (see 
e.g. jMota fc van de Bruckl2'ÖÖilMota fc Barrowl20Ö4allbl) l. 
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Figure 11. Evolution of 5G/G as a function of log(l + z) for 
overdensiUes which virialise at z = O in a Kcdm model. Thick 
solid line 2t^ + 3 = 4 X 10®, thin solid line 2üj + 3 = B'^cj)^'^P+^^ , 



dashed-line 2c^ + 3 = 2A\1 - ^l" 
C\ ln(^)|-4^ dotted line 2lj + 3 = 
is normalised in order to have Go 
densities. 



dash-dotted line 21x1 + 3 = 

1- {-Ê-Y\'^ ■ Each model 
Gc{z = 0) inside the over- 



Figure 13. Values ofG/G as a function oflog{l + Zy) at virial- 
isation. Upper plot: evolution inside the overdensity. Lower plot: 
evolution in the background universe. Thick solid line 2aj + 3 = 
4 X 10®; thin solid line Ilü + 3 = B^(j); dashed-line 2a; + 3 = 
2A|1 - ^|-2; dash-dotted line 2lj + 3 = C|ln(^)|-4; dotted 

line 2ti; + 3 = D|1 — (-f^)^l^-'-. Each model is normalised to have 
Go = Gc{z = 0) inside the overdensities . 
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Figure 12. Values of AG/G as a function of log(l + Zy) at 
virialisation. Upper plot: evolution inside the overdensity. Louier 
plot: evolution in the background universe. Thick solid line 2lü + 
3 = 4 X 10®; thin solid line 2üj + 3 = B^ipi dashed-line 2a; + 3 = 
2A\l-^\-'^: dash-dotted line 2a; + 3 = Cl ln(-^)|-''; dotted line 

2aj + 3 = D\l — (-^)^\~^ . Each model is normalised to have 
Go = Gc(z = 0) inside the overdensities. 



Difïerent initial conditions will lead to slightly different al- 
lowed values of u) and to different cosmological behaviours 
of G and Ge. 



6 CONCLUSIONS 

We have studied the inhomogeneous cosmological evolu- 
tion of the Newtonian gravitational 'constant' G within the 
framework of relativistic scalar-tensor theories of gravity, of 
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Figure 14. Values of SG/G as a function of log(l + Zy) at 
virialisation. Thick solid line 2a) + 3 = 4 X 10®; thin solid line 
2lü + 3 = B^if>; dashed-line 2üj -\-3 = 2A\1- 

line 2a; + 3 = C|ln(^)|-*; dotted line 2lü + 3 
Each model is normalised to have Go 
overdensities. 
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dash-dotted 

^2|-l 



ö|l-(^)1- 
Gc(z = 0) inside the 



which Brans-Dicke theory is the simplest and best known 
case. We began by first exploiting the conformal equivalence 
between these theories and general relativity to transform 
an existing solution of Einstein's field equations to a new 
exact spherically symmetrie inhomogeneous vacuüm cosmo- 
logical solution of the Brans-Dicke field equations. We then 
present a second spherically symmetrie perfect-fluid solution 
to the Brans-Dicke field equations which corresponds to an 
overdense inhomogeneity in an asymptotically homogenous 
and isotropic expanding cosmological background. These ex- 
act Solutions have simple enough form to allow the G{r, i) 
evolution to be investigated directly and can be used to 
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Figure 15. Values of SG/G as a function of Ac for a Acdm 
model. Thick-solid line 2a; + 3 = 4 X 10®/ solid line 2a; + 3 = -B^0; 
dashed-line 2üj + 3 = 2A|1 — -$-\~'^ : dash-dotted line 2üj + 3 = 

C|ln{^)|-4; dotted line 2c<; + 3 = D|l - {-^^f\~^- Each model 
is normalised to have Go = Gc{z = 0) inside the overdensities. 



will produce spatial inhomogeneities in G with a value which 
depends on the scalar-tensor theory used. While some mod- 
els produce higher values of G inside the overdensity, oth- 
ers produce a lower one. In spite of these differences, such 
spatial inhomogeneities are small SG/G ~ 10^''. The differ- 
ences in the time variations of G were found, with typically 
Gc/Gc ^ 10"^°s-^ inside the clusters and G/G < 10"^^s"^ 
in the background universe in the range of theories studied. 
Such variations in G do not significantly alter the virialisa- 
tion radius and the density contrast Ac of the virialised clus- 
ters from those in Standard Acdm and cdm models with con- 
stant G. Nevertheless, different scalar-tensor theories lead to 
slightly different results. There is also a dependence on the 
initial conditions. Different initial conditions will lead to a 
different value of ui and to different cosmological behaviours 
of G and Ge- Taken as a whole these analyses show that lo- 
cal observational limits on varying G made within our solar 
system or Galaxy must be used with caution when placing 
constraints upon the allowed cosmological variation of G on 
extragalactic scales or in the early universe. The universe 
is not spatially homogeneous and, in cosmological models 
where it can vary, nor is G. 



model the presence of a Schwarzschild-like mass in an ex- 
panding Brans-Dicke universe that approaches an idealised 
FRW universe asymptotically as r ^ oo. It is noted that 
far from the mass G varies very slowly with r and that as 
r — » oo the variation of G with r is removed altogether and 
only the cosmological evolution of G with time is relevant. 
Close to the mass, the r variation of G becomes more sig- 
nificant and we see that G ^ O at small r. We also note 
that in the limit tj ^ oo all space and time variation in 
G is removed from these solutions and general relativity is 
recovered. 

Next, we increased the complexity of the model by 
matching together two vacuüm FRW-Brans-Dicke universes 
of zero and positive curvature on a spacelike time slice. This 
results in a simple model for a spherically symmetrie pertur- 
bation in the density and curvature that is foliowed into the 
non-linear regime. The different evolutions of G in the two 
regions were determined and a comparison was made of the 
different evolutions of G on spacelike slices of differing time. 
We see from this toy model that the value of G does indeed 
have a different value in regions that have decoupled from 
the expanding cosmological background and coUapsed, com- 
pared to its value in the background itself. We were then able 
to repeat this construction for a more realistic matching of 
two FRW Acdm universes of different curvature in an arbi- 
trary scalar-tensor gravity theory. We foliowed the develop- 
ment of spherical overdensities through their expansion, sep- 
aration from the background, turnaround, and subsequent 
coUapse. Applying a simple approximation to virialise the 
coUapsing overdensities, we were then able to study the dif- 
ferences in G{x, t) between the non-expanding overdensity 
and the expanding background universe. We highlight as a 
special case the simple Brans-Dicke theory but also present 
results for other scalar-tensor gravity theories, specified by 
their defining coupling function a;(0). Although each theory 
predicts a different detailed cosmological evolution of G, a 
feature common to all of them is the sharp difference be- 
tween the value and time evolution of G inside bound over- 
densities and in the background universe. These differences 
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APPENDIX A: INHOMOGENEOUS FIELD 
EQUATIONS 

Substituting 1231 1 into Q gives, for a perfect fluid 



a f ^\ ^ f ^\ r,^ 't' 4> 
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as the r", T^^ and r°" equations, respectively. The 
propagation equation ij^J now becomes 



where fc = ^ / |t|^ ■ For e" and e^ of this form, equations 
lAll l. <A3l l and 1A4II become 
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and the only other non-trivial field equation is that for 



yi 



a cj> (f) (f> a (f) 2 (f) 

We now assume (j) is of the form (ji{r, t) = (j>(r)<j>{t) and 
look for Solutions to the set of equations 



Stt (3a;pe'' + 3pe'' + upe") 



mt) 



(2a; + 3) 



a a(t>{t) u) I (j){t) , 



Where 1A15I I was obtained by substituting IIAH and 
ljA3|l into l|A4fl and discarding the terms involving r deriva- 
tives, as these now sum to zero. We see that 1A13II . <A14t 
and lIAlSII are simply Qi © a-nd ^ with fc = O, Pfrw ~ 
pe" and Pfrw ~ pe", where subscript frw denotes a quan- 
tity derived from the field equations using the FRW metric. 
We also have, from Tt" = O, that 
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Looking for solutions of the form a (x t^ and (j){t) oc f" 

(A6) gives, on substitution into IIA13L the relation y — 2 — 3x7 

(assuming an equation of state for the Universe of the form 

p = (7 — l)p). Using 1A13L IIA14II . and this relation then 

gives the solutions 
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Such solutions are given by (INariailllQfiflf) as 
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These are exactly the same as would be expected for the 
scale factor and the BD field in a flat FRW Universe ([Nariail 
Il968h . The form of p is then given by (IA16t and (lAlOU as 



^('■'*)=^"(5)) 

and <l){r, t) is given as 
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